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INTRODUCTION 

It  has  been  suggested  that  the  direct  use  of  solar  energy  to  propel  spacecraft  from  low  earth 
orbit  to  geosynchronous  orbit  is  best  performed  by  absorbing  radiant  solar  energy  directly  into  the 
rocket  propellant  M.  The  absorbed  radiant  energy  would  be  converted  through  collisions  into 
thermal  energy,  heating  the  propellant  to  provide  a  high  exhaust  velocity.  The  solar  plasma 
absorber  has  the  advantage  that  it  would  be  able  to  reach  higher  propellant  temperatures  than  the 
alternately  proposed  blackbody  absorber.  Another  advantage  is  the  tendency  of  the  plasma  to 
trap  reradiated  energy,  thus  reducing  loss  and  increasing  efficiency. 

The  preferred  gas  for  spacecraft  propulsion  is  hydrogen,  since  as  the  lightest  element,  it 
would  provide  the  highest  specific  impulse.  However,  hydrogen  is  a  poor  absorber  of  solar  energy 
because  it  is  transparent  in  the  visible  and  near  infrared  region  of  the  spectrum  where  the  solar 
flux  is  greatest.  To  better  absorb  solar  energy,  the  hydrogen  propellant  would  have  to  to  seeded 
with  lightweight  absorbing  elements.  The  alkali  metals  would  be  excellent  candidates  for  seedants 
since  it  is  already  known  that  they  absorb  well  in  the  visible  and  near  infrared  regions  of  the 
spectrum  t2l.  Alkali  hydrides  would  also  form  and  contribute  to  absorption  PI.  In  fact,  it  is  likely  that 
previously  unobserved  species  MXHX  (where  M  is  an  alkali  metal)  may  also  form  and  contribute 
to  absorption3. 


EXPERIMENT 

In  order  to  study  candidate  absorbing  elements  and  their  hydrides  in  the  high 
temperature, high  pressure  conditions  that  would  exist  in  the  solar  plasma  thruster,  the  Plasma 
Spectroscopy  Cell  (PSC)  (Fig.  1)  has  been  constructed  at  OLAC  Phillips  Laboratory,  Edwards 
AFB,  California.  The  PSC  uses  concentric  gas  flows  and  resistive  elements  to  create  high 
temperature  (potentially  3000  K),  high  pressure  (potentially  100  atmospheres)  metal  vapors  which 
can  be  optically  studied  for  their  absorptions  properties. 

The  metals  are  loaded  in  a  crucible  surrounded  by  a  rhenium  resistive  heating  ring  in  a 
section  of  the  PSC  known  as  the  atomizer  or  evaporator  region.  A  neutral  entrainment  gas,  such 
as  helium  or  argon,  carries  the  metal  vapor  through  the  central  tube  A  into  the  region  known  as 
the  hot  zone.  The  hot  zone  is  heated  by  a  molybdenum  susceptor  which  surrounds  the  concentric 
tubes  A,  B,  and  shroud  tube.  Tube  B  carries  a  gas  which  can  be  chosen  as  neutral  or  a  reactive 
gas  such  as  hydrogen.  The  metal  vapor  jet  formed  in  the  central  region  of  the  hot  zone  as  it  exits 
Tubes  A  and  B  contains  the  molecules  of  interest.  Optical  access  to  this  area  is  available  through 
four  opposing  windows  on  two  perpendicular  axes. 

Absorption  experiments  carried  out  with  the  PSC  over  the  last  three  years  have  generally 
followed  the  schematic  shown  in  Figure  2.  Both  xenon  and  tungsten  lamps  have  been  used  as 
sources  to  pass  light  through  the  PSC  and  into  a  spectrometer,  In  several  experiments,  the  extra 
beam  created  by  the  beam  splitters  has  been  used  with  an  additional  spectrometer  to  obtain 
spectra  of  wider  spectroscopic  range  or  higher  resolution  in  a  region  of  particular  interest. 
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Figure  1. 

The  Plasma  Spectroscopy  Cell  [3] 
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Figure  2. 

Schematic  of  PSC  Experiment 


DATA 

Much  of  the  data  obtained  in  the  PSC  experiments  described  above  has  been  in  the  form  of 
images  384  by  578  pixels.  The  long  dimension  corresponds  to  wavelength,  the  range  of  which 
was  300  to  665  nm.  The  shorter  dimension  corresponds  to  spatial  position  across  the  width  of  the 
metal  vapor  jet.  Thus,  each  image  file  contains  578  cross  sections  of  intensity  vs  position,  and 
384  cross  sections  of  intensity  vs.  wavelength.  The  image  files  for  the  original  intensity,  l0,  of  the 
source  lamp  and  the  transmitted  intensity,  I,  through  the  metal  vapor  jet  were  used  to  calculate  the 
images  of  absorbance  according  to 


where  B0  and  B  are  the  appropriate  background  intensity  images. 

It  is  easier  to  reproduce  and  visualize  a  single  absorbance  intensity  vs.  wavelength  cross 
section  of  an  image  file  rather  than  the  whole  image  file.  An  example  of  such  a  cross  section  is 
shown  in  Figure  3,  which  shows  absorbance  by  sodium  at  1720  K  and  a  pressure  of  one 
atmosphere.  Prominent  atomic  lines  are  at  330  nm  for  the  3s-4p  transition  and  the  broad, 
saturated  sodium  D  line  at  589  nm  for  the  3s-3p  transition.  Centered  on  500  nm  is  absorbance 

by  the  B'n«  ~  transition  of  Na2.  Another  diatomic  sodium  absorbance  in  the  ridge  between 
550  and  560  nm,  believed  to  be  attributed  to  the  ~  lJ£.+  transition. 

Experiments  performed  with  a  mix  of  lithium  and  sodium  produced  results  such  as  the 
absorbance  plot  shown  in  Figure  4,  taken  at  a  temperature  of  1600  K  and  pressure  of  one 
atmosphere.  Again,  the  sodium  atomic  peaks  for  the  3s-4p  and  3s-3p  transitions  appear, 
whereas  the  atomic  2s-2p  transition  of  lithium  at  670  nm  is  off  the  right  edge  of  the  range  and  only 

appears  as  a  high  ridge  of  line  broadening  mixed  with  the  ~  transition  of  Li2.  The  ridge 

centered  on  500  nm  with  more  sharply  defined  bandheads  is  the  Bln.  -  X1!*  transition  of  Li2. 
Transitions  from  NaLi  and  Li  triplets  may  be  mixed  in  with  the  Li  and  Na  features  identified  so  far. 
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Figure  3. 

Absorbance  of  Na  at  1720  K  and  1  Atmosphere 


Figure  4. 

Absorbance  of  Na,  Li  Mixture  at  1600  K  and  1  Atmosphere 

Adding  hydrogen  as  a  reactive  gas  allows  the  possible  formation  of  hydrides,  e.g.,  LiH  and 
NaH.  Figure  5  shows  a  predominantly  lithium  absorption  spectrum  taken  at  1900  K  and  one 
atmosphere.  Residual  sodium  can  still  be  seen  to  contribute  atomic  absorptions  at  330  and  589 
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nm.  The  molecular  transitions  of  B-X  and  A-X  from  Li2  are  seen  around  500  nm  and  above  61 0 
nnm,  respectively.  The  lithium  atomic  line  for  the  2p-3d  transition  can  be  seen  at  610  nm  to  the 

right  if  the  ridge  believed  to  belong  to  the  ~  l5^+  transition  of  Li.  The  remaining  absorption 
bands  from  330  to  430  nm  are  in  the  right  range  for  both  c'n - -  absorption  of  Li2  and 
Axl*  -  xlI+  absorption  of  LiH. 


Figure  5. 

Absorbance  of  Na,  Li,  H  Mixture  at  1900  K  and  1  Atmosphere 

ANALYSIS 

General  references  such  as  Huber  and  Herzberg’s  book  on  Molecular  Spectra  W  can  be  used 
to  give  an  indication  where  molecular  transitions  should  occur  in  the  electromagnetic  spectrum. 
However,  the  absorption  spectra  obtained  in  the  PSC  experiments  may  have  overlapping  features 
contributed  by  various  transitions  from  the  same  and  different  molecules  in  the  chamber.  In  this 
mix  of  transitional  spectra,  it  is  more  difficult  to  identify  which  species  are  contributing  significantly 
to  overall  absorption  of  radiant  energy.  A  more  detailed  analysis  is  necessary  to  identify  the 
contributions  potentially  useful  for  the  absorption  of  solar  energy  in  the  solar  plasma  thruster. 

To  more  carefully  identify  absorption  spectra  from  specific  transitions  of  a  diatomic  molecule, 
certain  standard  techniques  of  computation,  based  on  molecular  theory,  have  been  developed. 
The  method  begins  with  the  potential  energy  curves  of  the  transitions  in  question.  The  potential 
curves  are  used  to  generate  accurate  eigenvalues  of  the  radical  Schrodinger  equation.  Both  the 
potential  curves  and  the  eigenvalues  can  then  be  used  to  calculate  the  intensities  for  ail 
rotational-vibrational  energy  levels  of  the  two  electronic  states. 

To  illustrate  this,  consider  the  triplet  transitions  of  Li2  which  should  appear  in  the  spectra 
range  of  the  PSC  absorption  data  {330  -  665  nm).  Figure  6  shows  the  potential  curves  for 

the  i%+.  1 3nv  2}X* ,and  23/7f  states  of  Li2.  The  lowest  state,  l3^.+>  is  an  RKR  potential  curve 
generated  by  Zemke  and  Stwalley  [5],  which  has  only  eleven  bound  states. 
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Figure  6. 

Potential  Curves  of  Selected  Li  Triplets  [5,6,7] 

The  1 37T,  and  t+  states  are  from  ab  initio  calculations  of  Schmidt-Mink,  Muller  and  Meyer  [6]. 
The  1 state  is  purely  repulsive  and  has  no  bound  states.  The  state  was  found  to  have 

54  bound  states.  The  portion  of  the  23/7,  state  shown  in  Figure  6  is  an  RKR  potential  curve 
generated  from  molecular  constants  provided  by  Xie  and  Field  [7].  It  will  be  necessary  to  merge 
the  RKR  curve  with  theoretical  long  range  values  before  eigenvalues  can  be  calculated. 

Often  when  the  intensity  of  molecular  transitions  are  computed,  the  Franck-Condon  Factors 
(FCFs)  are  calculated  with  the  assumption  that  the  electronic  transition  dipole  moment  can  be 
treated  as  a  constant  over  the  intermolecular  distance.  For  this  analysis,  the  full  functions  of  the 
electronic  transition  dipole  moments  were  used  to  calculate  the  Franck-Condon  Factors,  as  they 
depart  significantly  from  a  constant  at  short  range.  Figures  7  and  8  show  the  functions  used  for 

the  i3£„+  -  l }ng  and  l%+  -  2 3Z*  transitions  of  Li2l  respectively,  as  calculated  by  Ratcliff,  Fish, 
and  Konowalow  [8]. 

The  intensities  of  the  i3£«+  -  2 3Z}  transition  of  Li2  are  represented  by  the  Franck-Condon 
factors  plotted  in  Figure  9.  Intensities  were  only  calculated  for  the  rotationless  transitions  from  the 

11  vibrational  states  of  l3^  to  the  54  vibrational  states  of  2 3Zf  .  Further  calculations  will  be 
performed  to  include  rotational  energies  in  the  transitions.  The  intensities  can  then  be  binned  into 
small  wavelength  ranges  corresponding  to  the  578  pixels  in  the  wavelength  direction  of  the  PSC 
image  data  for  comparison  with  the  data. 

For  bound-free  absorption,  the  intensities  are  calculated  directly  in  atomic  units.  Figure  1 0 
shows  the  computed  intensities  for  the  l3^«+  -  i3^*  transition  at  a  temperature  of  2000  K. 
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Figure  7. 

Electronic  Transition  Dipole  Moment  for  i3^.  -  l *nt  0f  Li2  [8] 


Figure  8. 

Electronic  Transition  Dipole  Moment  for  i3^.  -  2 3z;  0f  U2  [8] 
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Figure  9. 

Rotationless  Intensities  of  the  i3^«+  -  Transition  of  U2 


FigurelO. 

Intensities  of  the  U2  131  *  ~  l3/7»  Transition  at  a  Temperature  of  2000  K 

The  jet  profiles  of  absorbance  vs.  position  contained  in  the  PSC  image  data  are  not  true 
indicators  of  the  amount  of  absorbance  per  unit  length  occurring  in  the  metal  vapor  jet.  To  obtain 
the  radical  profile,  it  is  necessary  to  perform  an  Abel  inversion  on  the  lateral  profile.  For  a  measured 
lateral  intensity  profile  N(x)  of  a  specific  wavelength,  the  inverted  profile  f(r)  is  given  by 
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«')  -  -  h 


N'(x)dx 
(x2  -  r2)1/2 

where  R  is  the  edge  of  the  lateral  profile  [9],  However,  first  performing  a  derivative  on  the  data 
tends  to  amplify  the  effects  of  noise.  A  method  recommended  by  Freeman  and  Katz  [10]  first 
rewrites  the  integral  with  the  derivative  on  the  outside  as 


fir) 


JL  £  [  N(x)xdx 
™  dr  J  (x2  -  r2)1' 2 


Next,  the  data  is  fit  to  a  curve  of  the  form 


m  =  J  Ctf1  ~  **)' 

1-0 

which  can  be  integrated  by  parts  to  give  the  Abel  inversion 

fir)  =  X  Ci  ( * 2  ~  '•2)i_1/2 

i 

where 

i  _  2 W 
‘  n(2 i)l 

Using  the  method  above,  the  data  in  Figure  11  was  fit  to  a  curve  of  eight  order.  Original  lateral 
data  profile  represented  by  small  circles  was  fit  to  lateral  data  profile  as  solid  line.  Dashed  line  is 
Abel  inversion  multiplied  by  100.  The  resulting  Abel  inversion  has  been  scaled  by  a  factor  of  100 
to  show  the  narrower  shape  of  the  radical  profile  as  compared  with  the  original  data. 


Figurell. 
Abel  Inversion 
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CONCLUSIONS 

High  temperature  lithium  and  sodium  experiments  in  the  Plasma  Spectroscopy  Cell  (PSC) 
have  generated  a  large  volume  of  data  which  can  be  analyzed  with  modern  computational 
techniques.  Calculations  to  date  have  shown  that  these  techniques  can  be  used  to  generate 
intensity  profiles  that  will  help  to  identify  the  specific  molecular  transitions  observed  in  the  PSC 
experimental  data  and  thus  help  to  characterize  the  mix  of  metal  vapors  to  be  used  for  solar  plasma 
propulsion.  Further  development  of  these  computational  methods  will  enable  the  direct 
comparison  of  the  calculated  absorption  profiles  to  the  PSC  data. 


REFERENCES 

1.  Larson,  C.W.,  Solar  Plasma  Propulsion.  AFRPL-TR-85-096,  Air  Force  Rocket 
Propulsion  Laboratory,  Edwards  AFB  CA,  1987. 

2.  Stwalley,  W.C.,  and  Koch,  M.E.,  “Alkali  Metal  Vapors:  Laser  Spectroscopy  and 
Applications",  Optical  Engineering,  Vol.  19,  No.  1,  Jan/Feb  1980,  pp.  71-84. 

3.  Larson,  C.W.,  The  Spectroscopy  of  Hvdroqen/Metal-Vapor  Mixtures  at  High 
Temperatures  and  Pressures.  AL-TR-88-080.  Astronautics  Laboratory,  Edwards  AFB 
CA,  1988. 

4.  Huber,  K.P.,  and  Herzberg,  G.,  Molecular  Spectra  and  Molecular  Structure  IV  Van 
Nostrand  Reinhold  Co.,  New  York,  1979. 

5.  Zemke,  W.T.,  and  Stwalley,  W.C.,  “Analysis  of  Long  Range  Dispersion  and  Exchange 
Interactions  of  Two  Li  Atoms",  Journal  of  Physical  Chemistry,  Vol.  97,  1993,  pp. 
2053-2058. 

6.  Schmidt-Mink,  I.,  MUIIer,  W.,  and  Meyer,  W.,  “Ground-  and  Excited-state  Properties 
of  Li2  and  Li2+  from  Ab  Initio  Calculations  with  Effective  Core  Polarization  Potentials”, 
Chemical  Physics,  Vol.  92,  1985,  pp.  263-285. 

7.  Xie,  X,  and  Field,  R.W.,  “Perturbation  Facilitated  Optical-Optical  Double  Resonance 
Spectroscopy  ofthe6Li233  2  g+,  23ng,  i3Ag,  b3nu,  and  a3  2  u+  States  ”,  Journal  of  Molecular 
Spectroscopy,  Vol.  117,  1986,  pp.  228-244. 

8.  Ratcliff,  L.B.,  Fish,  J.L.,  and  Konowalow,  D.D.,  “Electronic  Transition  Dipole  Moment 
Functions  for  Transitions  Among  the  Twenty-Six  Lowest-Lying  States  of  Li2”,  Journal  of 
Molecular  Spectroscopy,  Vol.  122,  1987,  pp.  293-312. 

9.  Tourin,  Richard  H.,  Spectroscopic  Gas  Temperature  Measurement.  Elsevier 
Publishing  Company,  New  York,  1966. 

10.  Freeman,  M.P.,  and  Katz,  S.,  “Determination  of  a  Radiance-Coefficient  Profile  from 
the  Observed  Asymmetric  Radiance  Distribution  of  an  Optically  Thin  Radiating  Medium”, 
Journal  of  the  Optical  Society  of  America,  Vol.  53,  1963,  pp.  1172-1179. 


10 


